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On  Recoverable  Internal  Energy  in  Linear  Viscoelasticity* 

by 

S.  Breuer**  and  E.  T.  Onat*** 

Summary.  A  linear  viscoelastic  solid  is  subjected  to  a  given 
deformation  history.  A  portion  of  the  work  done  by  the  stresses 
during  this  deformation  is  converted  into  heat,  while  the  remain¬ 
ing  portion  increases  the  internal  energy  (per  unit  volume)  of  the 
solid.  A  fraction  of  the  increase  in  internal  energy  can  be  re¬ 
covered  by  subjecting  the  solid  to  an  appropriate  future  deforma¬ 
tion.  The  paper  is  concerned  with  the  question  of  maximizing  the 
recoverable  energy  by  means  of  an  optimum  future  deformation. 

It  is  shown  that  the  determination  of  the  optimum  de¬ 
formation  requires  the  solution  of  an  integral  equation  of  the 
Wiener-Hopf  type.  This  equation  is  solved  in  the  case  where  the 
relaxation  modulus  is  given  as  a  sum  of  exponential  functions. 

The  maximum  recoverable  internal  energy  is  then  expressed  as  a 
functional  of  second  degree  of  the  given  deformation  history. 

It  is  observed  that  the  maximum  recoverable  energy  pro¬ 
vides  a  lower  bound  to  the  internal  energy  of  the  solid.  It  is 
hoped  that  use  could  be  made  of  the  concept  of  maximum  recoverable 
energy  in  studies  concerned  with  the  thermodynamics  of  linear  vis¬ 
coelasticity, 

*  The  results  communicated  in  this  paper  were  obtained  in  the 
course  of  research  sponsored  by  the  Office  of  Naval  Research 
under  Contract  Nonr  562(10)  with  Brown  University. 

**  Research  Assistant,  Division  of  Applied  Mathematics,  Brown 
University,  Providence,  R.  I. 

***  Professor  of  Engineering,  Brovn  University,  Providence,  R.I. 


t 


C 11-82 


2 


1.  Introduction.  Consider  Isothermal  deformations  of  a  linear 
viscoelastic  solid  in  simple  tension  or  compression.  Let  d(t) 
and  e(t)  denote,  respectively,  the  stress  and  infinitessimal 
strain  components  at  time  t.  We  adopt  the  stress-strain  relation 
in  the  form  [l]* 

o(t)  =  J  Q(t-r)  dx  ,  (1) 

-  00 

where  G(t)  is  the  relaxation  modulus,  which  need  be  defined  only 
for  non-negative  values  of  its  argument. 

Let  the  material  be  subjected  to  a  given  deformation 
e(t)  in  the  past,  i.e,  in  the  time  interval 

-  co  <  t  <  0  , 

where  t=0  refers  to  the  present  instant.  The  work  done  by  the 
stress  in  this  interval  is  given  by 

W  *  f  <j(t)  dT  ,  (2) 

-co  -co 

By  using  the  constitutive  equation  (1)  and  by  extending  the  range 
of  definition  of  G(t)  with  the  relation 

G( — t )  *  G(t)  ,  (3) 

(2)  can  be  written  in  the  following  form, 

*  Numbers  in  square  brackets  refer  to  the  bibliography  at  the 
end  of  the  paper. 
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0,00 

w  s  2  f  /  GCt-n:)  v(t)  v(t)  dt  dT  ,  (4) 

-oo  -00-00 


where  v(t)  denotes  the  rate  of  strain  at  time  t, 


V(t)=f.  (5) 


It  was  shown  in  [2]  that  the  second  order  functional 
o 

W  is  positive  definite  if  G(t)  is  a  steadily  decreasing  func- 
-co 

tion  of  time  (for  positive  times)  which  is  convex  from  below 

and  tends  to  a  non-negative  asymptotic  value  as  t  — >oo.*  In 

the  present  paper  we  shall  be  concerned  with  the  relaxation 

moduli  which  possess  the  above  mentioned  properties* 

Since  viscoelastic  solids  are  dissipative,  a  certain 
o 

portion  of  the  work  W  will  be  converted  into  heat  during  the 

-oo 

course  of  the  deformation,  while  the  remaining  portion  will 
increase  the  internal  energy  of  the  solid.  We  now  consider  the 
possibility  of  converting  the  internal  energy  of  the  solid  into 
useful  work  by  purely  mechanical  means.  For  this  purpose,  we 
subject  the  solid  to  further  deformation  in  the  future,  i.e,  in 
the  time  Interval 


0  <  t  <  co  , 


♦ 

Less  restrictive  conditions  on  G(t)  which  still  ensure  positive 

o 

definiteness  of  W  have  been  discussed  in  [3]. 

-00 
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GO 

and  denote  by  W  the  work  done  by  the  stress  In  this  Interval, 
o 

If  the  future  deformation  Is  chosen  in  such  a  way  that 


oo  oo 

W  ■  ^  d(t)v(T)  d t  <  0  ,  (6) 

then  it  is  said  that  useful  work  is  extracted  from  the  material 
in  the  interval  0  <  t  <  oo.  In  other  words,  deformations  in 
0  1  t  <co  satisfying  (6)  enable  one  to  recover  a  portion  of 
the  internal  energy  of  the  solid  which  has  been  subjected  to  a 
given  deformation  in  -oo<t  <0,  the  recovered  portion  of  the 
internal  energy,  Er,  being  given  by  the  obvious  relation 

no 

E  -  -  W  .  (7) 

o 


The  present  paper  is  concerned  with  the  following  prob¬ 
lem:  Given  v(t)  in  the  interval  -co<t  <  0,  what  function  v(t), 
defined  in  the  interval  0<t  <oo,  maximizes  the  recovered  in¬ 
ternal  energy  E  .  The  corresponding  maximum  value  of  the  re- 

o 

coverable  internal  energy  will  be  denoted  by  ErB). 

-oo 

o 

It  will  be  seen  that  E  will  depend  on  the  relax¬ 


ation  modulus  of  the  material  and  on  the  given  strain  history  in 

-oo  <t  <0.  If  one  assigns  the  value  zero  to  the  internal  energy 

o 

of  the  unstressed  virgin  material,  then  Erm  will  provide  a 


of  the  solid  at  the  termina¬ 


tion  of  the  application  of  the  given  deformation  in  -oo<t<0. 
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This  remark  follows  from  the  observation  that  the  optimum  ex¬ 
traction  process  in  0<t  <co  will  be  accompanied,  in  general, 
by  further  dissipation.  It  is  hoped  that  the  availability  of  a 
lower  bound  to  the  internal  energy  of  solids  may  be  useful  in 
the  study  of  thermodynamics  of  linear  viscoelastic  solids* 

In  section  2  it  is  shown  that  the  above  problem  of 
maximization  leads  to  an  integral  equation  of  the  Wiener-Hopf 
type  for  v(t)  in  0  <  t  <co.  Section  3  is  devoted  to  the  solution 
of  the  integral  equation  in  the  case  where  G(t)  is  given  as  a 
sum  of  exponentials.  It  will  be  seen  that  in  this  case,  the 
optimum  straining  which  maximizes  Ep  involves  a  sudden  applica¬ 
tion  of  strain  at  t=0,  followed  by  strain  rates  which  decrease 

in  magnitude  exponentially  with  increasing  time.  Section  4  is 

o 

concerned  with  the  evaluation  of  Erm  and  comparison  of  it 

-co 

with  the  free  energy  of  the  solid. 

It  may  be  appropriate  to  note  here  that  although  the 
discussion  in  the  present  paper  is  restricted  to  the  case  of 
simple  tension  or  compression,  there  is  no  essential  difficulty 
in  extending  the  present  concepts  and  analysis  to  the  general 
deformations  of  a  viscoelastic  solid. 

2.  Maximum  recoverable  internal  energy.  Basic  integral 
equation.  Using  the  definition  of  work,  the  recovered 
internal  energy,  Ep,  can  be  expressed  in  the  following  manner, 


3 
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6 


E 


r 


oo  o  oo 

w  =  w  -  w  , 

o  -oo  -oo 


o 

where,  according  to  (4),  W  depends  entirely  upon  the  given 

-oo 

deformation  in  -oo  <t  <0,  but 

W  *  i  f  J  G(t-'c)v(t)v('c)dtdT  ,  (9) 

-oo  -oo  -oo 


and  hence  Er,  depend  upon  the  entire  history  of  strain.  In  view 

of  (8)  the  problem  of  maximizing  Ep  reduces  to  the  following 

problem;  Given  v(t)  in  the  interval  -co  <t  <0,  to  find  v(t) 

oo 

in  O^t  <oo  which  minimizes  the  total  work  W,  i.e.  (9)* 

-oo 

Let  v(t)  be  the  function  minimizing  (9).  Consider  the 
set  of  functions  given  by 


u(t)  =  v(t)  +  ecjCt )  ,  -oo  <t  <oo  ,  (10) 


where  e  is  a  real  parameter  and  w(t)  is  arbitrary  except  that 


u(t)  b  0  ,  for  -co  <t<0.  (11) 

CO 

Denoting  by  W  [u]  the  work  done  by  the  stress  during  the  time 
-co 

interval  (-oo,  co)  in  the  course  of  the  deformation  character¬ 
ized  by  u(t),  we  obtain  from  an  expression  similar  to  (9), 

oo  co  _  co  op  co 

W[u]  -  W[v]  =  W[u>]  +  eT  f  0(t-T)v(t)u(T)dt d*r  • 

-oo  -co  -co  -&  -db 


(12) 
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oo 

Since  W  [w]  is  positive  definite  for  the  class  of  relaxation 
-co 

oo 

moduli  considered  in  this  paper,  (12)  shows  that  W[v]  will  be 

-oo 


a  minimum  provided  that 


co  oo 

Jk  u(t)  £  0(t-T)v(T)dTdt  *  0  .  (13) 


On  the  other  hand,  since  w(t)  is  arbitrary  for  t£0  and  zero 
otherwise,  we  must  have,  if  (13)  is  to  hold, 


oo  • 

|  0(t-t)v(t)dT*  0  ,  0  <  t  <oo.  (14) 

-co 


Moreover,  since  v('c)  is  known  for  -co  <  t  <0,  (14)  reduces  to 
f(t)  -  |  G(t-T)v(T)dT  ,  0£t  <co  ,  (15) 


where  f(t)  is  a  known  function,  defined  by 


o 

f(t)  »  -  fG(t— r)v(T)dT  ,  0<(t  <co.  (16) 


Equation  (15)  is  the  desired  integral  equation  for  the 
unknown  function  v(t)  in  0<t  <co  which  maximizes  E ^ ,  the 
_  recovered  internal  energy.  (2.5)  is  of  the  Wienpr-IJopf  type  and 
it^  solution  r>ay  be  arrived  at  by  the  Wiener-Hopf  technique 
[4],  In  the  present  study,  however,  we  shall  deal  with  a  partic¬ 
ular  class  of  relaxation  moduli  for  which  the  solution  of  (15) 
may  be  obtained  by  more  elementary  means. 
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Before  going  on  to  the  discussion  of  this  case,  we 

o 


shall  establish  a  result  concerning  Erm.  For  this  purpose, 

-oo 


let  us,  by  using  (3),  rewrite  (8)  in  the  following  form, 


E  s  -  fv(t)  J G(t-T)v(T)dTdt  -*  J*  f  G(t-T)v(t)v('c)dtd,&. 
o  -co  o  b 


(17) 


If  we  now  regard  v(t),  0<:t  <co,  as  the  solution  of  (15),  then 

o 


the  left  hand  side  of  (17)  becomes  E  ,  while  the  first  term 

-co 


on  the  right  hand  side,  in  view  of  (14),  takes  the  form, 


UJ  uu 

7  v(t)  I 


G(t-T)v(T)dTdt  , 


so  that  (17)  reduces  to 


=  1 
rm  “  2 


00  00 


n 


G(t-T)v(t)v(T)dtdT  . 


(18) 


Equation  (18)  shows  that  Erm  is  a  positive  definit( 

-oo 


functional  of  the  optimum  extraction  process.  This  implies 


that  for  the  class  of  materials  considered  in  this  paper  it  is 


always  possible  to  extract  useful  work  from  a  body  which  is 


subjected  to  non-trivial  (i.e.  non-zero)  deformations  in  the 


past. 


We  also  note,  for  future  reference, that  (18)  together 
with  (15)  yields 
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2  EM  =  J  f(t)v(t)dt  ,  (19) 

-00  0 

where  v(t),  0£t  <oo,  is  the  solution  of  (15). 


3,  Solution  of  the  Integral  equation.  We  shall  restrict  our 
attention  in  this  section  to  the  class  of  relaxation  moduli  de¬ 
fined  by 

N  -a, |t | 

G(t)  *  2c,e  1  ,  (20) 

i=l  1 


where  c^  and  a^^  are  positive  constants  and  a^a^.  For  con¬ 
venience,  we  define 

o  ait 

Kj  -  J  e  1  v(t)dt  ,  i*l,2,...,N  ,  (21) 

-oo 

where  v(t)  is  the  given  strain  rate  history  defined  in  -oo  <t  <0. 
By  combining  (20)  and  (21)  with  (16) ,  the  basic  integral  equation 
(15)  can  be  written  in  the  following  form, 

N  -a,t  oo 

-  2  c^e  «  |G(t-T)v(n:)dT  ,  O^t  <co, 

(22) 


where  G(t)  is  given  by  (20), 

In  seeking  a  solution  of  (22)  we  may  adopt  the  follow¬ 
ing  reasoning.  Suppose  N=1  in  (20),  Then,  as  it  is  well  known, 
the  viscoelastic  material  is  a  Maxwell  body,  represented  by  a 
combination  in  series  of  a  spring  and  a  dashpot.  The  internal 
energy  of  the  material  is  stored  in  the  spring,  and  it  can  be  re¬ 
covered  by  a  sudden  displacement  which  restores  the  spring  to  its 
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original  length.  Mathematically,  this  observation  implies  that 
the  solution  of  (22),  in  the  case  of  N=l,  is  of  the  form 

v(t)  *  A  6(t),  (23) 


where  6(t)  is  the  Dirac  delta  function*  and  A  is  a  constant  re¬ 
lated  to  the  past  deformation  of  the  body.  On  the  other  hand, 
in  the  case  of  N=2,  (22)  becomes 

-a-,t  -a«t  t  -a,  (t-T)  -a5(t-T) 

-c^K^e  -c^e  =  J  [c^e  +c2e  d  ]  v(,c)dT 


00  a,  (t-T)  a9(t-T) 

+  J  Lc-je  +c2e  c  ]  v(T)dr  . 


(24) 


The  integrals  on  the  right  hand  side  of  (24)  may  be  eliminated 
by  repeated  differentiation  of  (24)  with  respect  to  t.  Carrying 
out  this  procedure  one  obtains  a  differential  equation  for  v(t) 
which  in  turn  leads  to  the  solution  of  (24)  of  the  form 

v(t)  =  B5(t)  +  Ce"Xt  ,  (25) 

where  B,  C  and  X  are  suitable  constants. 

The  form  of  (25)  now  suggests  that  in  the  general  case 
of  (22)  we  may  try  a  solution  of  the  form 

N  -b,t 

v(t)  =  knb(t)  +  2  A,e  2  ,  (26) 

1  j*2  2 

For  the  present  purposes  the  Dirac  delta  function  is  defined 
in  the  following  manners 

00 

6(t)  *  0,  vfcen  t^O,  and  J  5(t)dt  *  1. 

0 
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where  the  constants  ^ ,  b^  are  to  be  determined  by  further 
considerations.  For  this  purpose,  we  introduce  (26)  into  (22) 
and  obtain 


N  -a.t  N  -a.t  N 

2  c^K.e  =  An  2  c.e  -  2 
1*1  1  1  1  i-1  1  j-2 


+  2 


N 

2  A.e 
j=2  J 


-V 


c^e 


-a^t 


N 

i.i 


N  cia1 
1=1  al'bJ 


(27) 


where  b,  are  subject  to  the  obvious  requirements 


b,  >0  ,  b,  4  ^  .  (28) 


We  now  demand  that 

N  c.a. 

s  -J-V  “  0  ,  3-2,3,...,  H,  (29) 

i-i  a^-bj 


2 

i.e.  we  desire  to  choose  bj  as  the  N-l  roots  of  the  equation 

N  Cja, 

2  4-±-  *  0.  (30) 

i=l  a^-  x 

p 

It  is  clear  that  b*  determined  from  (30)  will  satisfy 


»3  4  \  •  (31) 

Moreover,  it  can  easily  be  shown  that  (30)  possesses  N-l  distinct, 
positive  roots,  so  that  (28)  will  be  fulfilled  if  b^  are  taken 
as 

bj  *  +  >/xj  ,  Je2,3».»»>  N,  (32) 


I 
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where  x.,  are  the  roots  of  (30), 

J  -a.t 

Using  (29)  and  observing  that  e  are  linearly  in¬ 
dependent,  we  obtain  from  (27)  the  following  system  of  equations 


for  Aa, 


■Ki  =  ai'J2vV  111,2 . 


(33) 


It  is  shown  in  the  appendix  that  the  solution  of  (33) 


is  given  by 


n  (ai‘V 

m-n  A  8 


A1  -  -  2  Ki  ¥ - 


w 


h * 


n  tb -b  )  1=1 

q=2  q  3 


,!,Ki  N 


N  N 

n  (a  -bJ  n  (a.-b_) 
p=l  p  3  s*2  1  3 

_ 


,  ja2,3,* ,  N( 


n  (a.  -  a_) 
r=l  1 


(35) 


Equation  (26),  supplemented  by  (29),  (32),  (3*0  and 
(35),  determines  the  solution  of  (22)  explicitly. 

Equation  (26),  in  conjunction  with  (5),  reveals  that 
the  optimum  loading  program  consists  of  a  sudden  displacement 
followed,  for  N>1,  by  the  straining  program  characterized  by  the 
sum  of  exponentials  in  (26). 
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4.  Evaluation  of  the  maximum  recoverable  Internal  energy. 

Comparison  with  free  energy.  We  may  start  the  evaluation  of 

o  . . . . . ~  ’ 

E  rm  ty  combining  (19)  and  (26)  with  (16),  (20)  and  (21)* 

-oo 


o  oo  N  -b,t  N  -a<t 

2  B  ■  -  f  [A,6(t)+  E  U  3  ]  2  c-K-e  1  dt.  (36) 
-oorm  0  j*2  3  1=1  1  1 


Carrying  out  the  integration  in  (36)  we  obtain 


0 

2  E. 


N 


N 


N  A 


-A-i  2  CjKj  -  2  CjKj  2  •"  jv 
1  1=1  1  1  1=1  1  1  J=2  ai+b 


-00 


rm 


’J 


(37) 


Now  we  make  use  of  (34)  and  (35)  to  express  (37)  in  terms  of  K^, 
After  some  manipulations,  given  in  the  appendix,  we  obtain 


where 


and 


0  N  N 

E  ■  2  2 

-co™  i=l  j 


=1  ’ 


N  ? 

n  (a4+a_r 
r*=l  1  r 

_  rjfi _ 

"’ll  ~  N  N  « 

2  2  c  a  n  (at+b_)^ 

p»l  P  P  ss2  1  3 


aici 


>  0  • 


N 

A(ai+ar)(aj+ar) 

=  clc.1aia.1(al+a1)  r^l.,1 _ 


ij'ji  "  N 

2  c  a^ 
p-1  p  p 


N 


II  (a,+b  )(a4+b  ) 

.=0  1  B  Jo 


ss2 


(38) 


(39) 


>0,  i/j. 


(40) 
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o 

In  view  of  (21),  we  see  that  E  is  the  following  second  order 

-cor 

functional  of  the  given  strain  history  in  -oo<t<0, 

o  o  o  N  N  a.t  + 

Erm  =  f  J  !  2  V  3  v(t)v(>t)dt  it. 

-oo  -oo  -co  1=1  j=l  J 

(41) 


For  the  class  of  materials  considered  in  this  section, 
it  is  possible  to  determine  the  elastic  energy  stored  in  the 
material  by  regarding  the  solid  as  a  network  of  linear  elastic 
and  viscous  elements  (Staverman  and  Schwarzl  [5]  ,  Bland  [6] 

and  Hunter  [7]).  These  authors  deduced,  in  independent  works, 
the  following  expression  for  the  elastic  energy  stored  per  unit 
volume  of  the  material  subjected  to  simple  tension  or  compression,  . 


0  "|  0  0 

F  a  2  I  I  OCt+t)v(t)v(T)dt  dT  ,  (42) 

-oo  -oo  -do 

o 

where  F  denotes  the  elastic  energy  stored  at  the  termination 
-oo 

of  the  given  deformation  v(t)  applied  in  -oo<t<0,  and  where  G(t) 

is  given  by  (20).  Since  the  deformation  has  taken  place  under 

o 

constant  temperature,  it  is  appropriate 'to  interpret  F  ,  from 

-co 


the  point  of  view  of  thermodynamics,  as  the  free  energy  in  the 
sense  of  Helmholtz. 
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We  note  that  (42)  can  be  written  in  the  follow 
ing  form  by  making  use  of  (20)  and  (21), 


N 


(43) 


Since  the  extraction  of  internal  energy  discussed  in 
the  previous  sections  entails  dissipation  due  to  viscous  effects 
we  expect  the  following  inequality  to  hold, 


(44) 


Indeed,  as  shown  in  the  appendix,  by  using  (38),  (39),  (40)  and 
(43),  we  obtain 


0 

F 

-00 


N  N  9 

2  2  cl  ^(Kj-KJ^  . 

1=1  j.i  1  3 


(45) 


Since  >  0  [cf.  (40)],  (45)  shows  tfctt  (44)  is  fulfilled. 
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Appendix 

Our  first  task  is  to  solve  the  system  (33),  which  we 
repeat  here  for  convenience 

N  A, 


-  K,  = 


i  s  Ai  “  2  a  :b  ,  1=1, 2,...,  N. 
1  1  j=2  ai  Dj 


(46) 


Denote  the  determinant  of  the  system  (46)  by  A.  Then 

.  1 


A  -  (-1) 


N-l 


al"b2  V'fo 


a2‘b2 


a2_bi( 


Vb2 


•  •  • 


.  (47) 


Next,  subtract  the  first  row  of  the  determinant  in  (47)  from  each 


succeeding  row  to  obtain 


H 


where 


k,N 


N 


II  (a,-bj 

j=2  1  J 


2,N  , 


(48) 


1 

1 

1 

ak"^k 

Wi  "• 

ak-bN 

1 

1 

1 

^+1“^ 

wVi  - 

^+1"^ 

1 

1 

1 

“N'hk 

Wl  - 

2^c<N 


(49) 
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(49*) 


In  (49)  we  now  subtract  the  first  row  from  each  succeeding  row 
to  obtain 


N 


II  (a^-a  ) 
A  ^  r=k+l  K  r 


k,N  N 
n 


P=k 


1 

1 

1 

1 

1 

1 

^+1"^ 

^rVi 

ak+r% 

1 

i 

i  . 

ak+2_bk 

ak+2"bk+l 

1 

1 

1 

Wl 

V%  i 

(50) 


Finally,  in  the  determinant  in  (50)  we  subtract  the  first  column 
from  each  succeeding  column  to  obtain 


k,N 


N  N 

n  (ak-a  )  n  (b  -tO 

r°k+l  K  r  g.ktl  8  K  A 
N  N  k+l,N 


(51) 


n  (eu^-b  )  n  (a  -tO 


p=k  v  q=k+l 
Hence,  by  (48)  and  (51)  we  conclude 

N-l  N  N-l  N 

IT 

A  =  (■ 


in  -x  n  in -x  in 

n  n  (ai-aJ  n  n  (VV 

=  (.1)N-1  ^■llLi-L.Jgg-.gyl -1  p  4  o  . 


N  N 
n  n  (a  -b  ) 

r=l  s=2  r  3 


(52) 
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Next,  from  (46)  we  have 


K1 

1 

1 

1 

% 

arb3 

t  •  • 

K2 

1 

a2-b2 

1 

a2-b3 

Ml 

a2“^ 

% 

1 

1 

1 

®N"b 2 

£ 

1 

O' 

C«J 

•  •• 

(53) 


Now  in  the  determinant  In  (53),  the  coefficient  of  is  seen  to 
be  n  defined  In  (49),  while  the  coefficient  of  Kj,  for  J>1, 
is  equal  to  the  expression  obtained  from  -Ag  N  upon  replacing 


by  al*  Usin§  these  facts  we  obtain 

N  J5,  (al'V 


A,  g  -  Z  K 


sg2  A  ° 

i  N 

n  (a1-a  ) 
r-1  1  r 


(54) 


which  establishes  (34). 

Similarly,  (46)  gives, 


Cll-82 


1  al’b2  al‘bl-l  Kl  al'bi+l  al'bK 

A  =  t-l)”*1  ,  _1_  _J _  K  1  _1_ 

J  A  a2”b2  a2"bj-l  2  a2"bi+l  ”* 


1  Vb2  *"  Vhj-l  vte”,aN-^ 


js2j3»*»»  »N« 

(55) 

In  the  determinant  in  (55)*  the  coefficient  of  K^,  1*1,2,...,  N, 
is  given  by  («l)i+^  A1’^,  where  A*’^  is  equal  to  the  determinant 
obtained  from  (-1)N~^A  upon  striking  its  ith  row  and  jth  column, 
A  being  given  by  (47).  Using  this  observation  we  find,  after 


some  calculations, 


II  (a  -b.)  n  (a.-b  , 
p=l  P  J  s=2  1  s 


Aj  =  '  N 


S  K, 


,  J*2,3». . » ,N, 


",  <W 


n  (a4-a  ) 
r*l  1  r 


(56) 


Thus  (35)  is  established. 


In  order  to  establish  (38),  we  first  show  that  the 

2 

numbers  bj  defined  in  (29),  necessarily  satisfy  the  relations 


9  0 

c  .a,  Varv 

N  **  N  0  9 

V  _  -  TT  /_2 


2  c  a_  II  (al-aj) 
p-1  p  p  r=l  1  r 

r^i 


(57) 
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of  (58)  by  ^(a^-b^),  f0P  some  to  obtain 


N 


1^2(arbk) 

JL  °iai  .  isL i _ 

Di  a^-b?  »  2  2, 

1  J  II  (a, -a*) 

r=l  1  r 

i^l 


*  ,N$  j=2,3,. .  •  ,N. 


The  right  hand  member  of  (63)  expands  as  follows, 


(63) 


N  2  ? 

k/.l _ 


N 


N 

=  -  2 


N  2  ? 
n  (af-bj) 

k**2  8  K 
_ 


N 


"  i  c  _a 
s/1  3  1 


,2_.2\  S=1  H  /  2  2\ 


5 L”(aK) 


r=l 

r^i 


r=l  3  r 
r^s 


(64) 


the  last  equality  following  from  (63).  Finally,  (62)  end  the 
combination  of  (63)  and  (64)  yield  the  system  of  equations 


I  1  °lal 


l=l°ll^  =  °  ’  3=2)3 . N- 


S  ciai  , 

2  r\  =  1  , 

1*1  Di 


(65) 


containing  N  equations  for  the  N  unknowns  gk  Denote  the  deter 
minant  of  the  system  (65)  by  A1#  Then 
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I 


1 

arb2 

1 

Xh5  ••• 

a2"b2 

1 

fi2  v2 

Vb2 

a'  -  pivp 

1 

al"b3 

1 

.2  h2 

a2  b3 

1 

.2  v2 

Vb3 

1 

77 

1 

77  ••• 

a2-bN 

1 

“n  ’to 

1 

1 

1 

We  see  from  (66)  that  A*  is  proportional  to  A,  defined  in  (47), 
provided  we  replace  a^  and'b^  in  (47)  by  aj  and  b y  respectively. 
Hence  AVO,  since  A^O  by  (52),  and  (65)  has  a  unique  solution.  Now 
from  (29)  we  see  that  the  first  N-l  equations  in  (65)  will  be 
met  if  we  put 

D1  a  D2  a  *•••  88  %  s  D  »  (67) 
where  D  is  any  constant.  On  the  other  hand,  the  last  of  (65) 
will  be  satisfied  if  we  put 

N 

D  e  2  cn  an  *  (68) 

p°l  *  v 

This  completes  the  proof  of  (57)# 

We  proceed  to  establish  (38),  It  is  obvious  from  (34), 
0 

(35)  and  (37),  that  Erm  may  be  written  in  the  form 

•00 
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o 

E 

-CO 


nn 


N  N 
2  2  a 

1=1  J«1 


^ulj  Ki  Kj  ’ 


(69) 


where  the  coefficients  a^j  must  be  determined. 

N  Aj 

Consider  the  expression  2  T-rzr  ,  for  some  k.  By 

j=2  V  D;j 


(56)  we  have 


N 


N  N  D  (a  -bj 

w  n  (a.-b  )  «  p=l  p  " 

2  JLL-=  -  £  K  g3  l  5l 

5=2  VbJ  WJ‘S  2  - 


N  A. 


r=l 

r^i 


Equation  (70)  can  be  rewritten  as 


jyvv  J=2  Wj^VV 

q^5 


(70) 


N 


N 

II  (aj-aJ 


N  A 

2  r- 

j=2  ak 


_  .  11  vo«rai, 

II  (a*-b  )  p=l  P  k 

1  S  1/ . 


1  N  "  ^ai‘V 

k-  -  ?  ^ - [ 

«  *-*  n  _  \ 


K 

n,  ^ai"ar^  n  W+ba) 
r=l  1  r  s=2  K  s 

iVi 


-  1]  , 


<71; 

which  may  be  verified  by  expanding  the  term  contained  in  square 

brackets  in  (71)  in  terms  of  .  From  (71)  we  find  that 

*  N 

N  ,  n  (a+aj 

A  v  D  (ai"^Q^  P=1  P 

--*•*!* -LLtf - x] 


N 


N 


2  c^Kj  2  —ir 
1=1  1  1  5=2  ai+b5 


1=1 


H  (a,-a  )  n  (a,+b J 

1  r  _ o  1  S 


r*l 

n*i 


s=2 


N 


N  N 

-  2  CjKj  SKj 

i*l  1  1  J-l 
^1 


N  S  (a  +a* ) 

n  (a.-b  )  p=l  P  1 

- 1]  .  (72) 


5  N  1  N 

Q ,(*l"ar)  n  (&4+b  ) 
r*l  3  r  s*2  1  s 
r^j 
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Similarly,  we  find  from  (54)  that 


A*i  2  CjKj  ~ 
1  i=l  1  1 


N  n  (ai"V 

**  y  c  k2  ' 
AciKi  N 

1‘1  H  (araj 

r=l  1  r 

n/i 


n  n  “  (arV 

-  &eiKi  *  k5  5? - 

3^i  r=1(a3‘ar) 

P?J 


(73) 

Combining  (72)  and  (73)  and  taking  account  of  (37),  we  reach, 


N  n  (a  +at) 

.  „  n  (8,-b  )  p=l  P  1 

2  E  =  Z  c,K?  ■8jj^ - P& - 

rm  ,  ,  i  i  N  N 

'®  1-1  n  (8,-aJ  n  (a.+b  ) 

r»l  1  r  <f2  3 

r^i 

N 

N  II  (a  +a.) 

N  N  “  (aJ-b s>  Z] 

♦,!,  ^ — 

#  £(arar>  “2{ai+v 


(74) 
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At  this  point  we  make  use  of  the  relations  (57)  satis-* 
fied  by  a^,  bj  and  c^.  If  we  employ  (57),  we  are  able  to  write 
(7*0  in  the  alternative  form 


o 

2  E 


N 


1  £  c?a,K?  ^ 


H 

II  (a.+a_)‘ 
r=l  1  r 


„rm  "  N  wiai“l  N 

■°°  Sea  11 
p=l  p  p 


n  (a.+b  )d 
s*2  1  3 


JL 


N  N 


N 
2 

p=l  P  p 


N 

ini(a1+ar)(aJ+ar) 

SGU _ 


A  Acicjaia3(ai+aj)KiKJ  N 

Vi  syw<w 


(75) 

If  we  now  define  and  cc^  as  in  (39)  and  (40), 
respectively,  we  see  that  (75)  becomes  identical  with  (38). 

Thus  (38)  is  established. 

It  remains  to  establish  (45).  To  this  end,  we  observe 
from  (38)  and  (43)  that 


0  0  N  5  N  N 

2F-  2E  =  2(0.-2044)^-22  2  o.iK.K.. 

-co  -oor  i=l  1  11  1  i=l  J*1  13  1  3 

3/1 


By  (39)  we  have,  in  view  of  (57), 


(76) 
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ci“2aii 


n  (a^+a  )  r  N 
P=1  1  p  n 

P^i _  s= 


n  (a4+b  ) 
.-0  *  3 


'1  N 


n  (a,+b  )  n  (a,+a  )  n  (a, -a 

S=:P  1  3  n=l  1  P  1  P 


n  (a1-b  ) 

f2  1  8 

N 

n  (84—  a_) 


Expanding  the  bracket  in  (77),  we  obtain 


(77) 


n  (a,+a  ) r  N 

^  K  11  (Va.1) 

Cl'2aii=  °i  ^ -  1  +  -E2~H - 

S22(al+V  Vl  (al+aJ>r21Ur-aJ) 

r/i,j 


II  (b  -aJ 

3*2  ,  3 

N 


N  JJ2(afV 

-1-2  - ^ - 

j=l  N 

j/i  (a.-aj  II  (a,-a  ) 
1  j  r-i  j  r 

r^i  J 


(78) 


Combining  the  terms  in  the  brackets  in  (78),  we  get 


cr2aii 


II  (a,+a  )  H 

P=!  P  N  n  (a  -b  ) 

2c  a  _  v  _L_  s°2  j  8 

n  (ai+V  &i  3  n  (a.-aj 

s*2  1  s  r=l  J  r 


(79) 
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By  (57)  and  (40),  (79)  yields 


cr2au  = 


Combining  (80)  and  (76)  we  finally  reach 


o  o 


2  F  -  2  E 
-oo  -oo 


rm 


=  2 


N 

2 

i=l 


(80) 


.  (81) 


since  by  0+0),  we  see  that  (81)  is  identical  with  0+5). 
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